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摘 要 : 中 亚 高 山 冰 川 区 地 形 复杂 ,站 点 观测 和 传统 实地 测量 范围 十 分 有 限 。 卫 星 激光 测 高 技术 可 实现 大 范围 冰 
川 表面 高 程 变化 监测 。 以 2003 一 2009 年 ICESat 激 光 测 高 数据 为 数据 源 ,参考 2000 年 的 SRTM 高 程 数据 ,建立 冰川 


区 点 云 去 品 及 其 精度 优化 算法 和 多 尺度 冰川 区 表面 高 程 时 空 变化 分 析 模 型 ,并 分 析 了 2003—2009 年 间 中 亚 山 区 整 
体 与 各 分 区 冰川 表面 高 程 时 序 变 化 。 结 果 表 明 :中 亚 高 山 冰 川 区 的 冰川 表面 平均 高 程 整体 呈 下 降 趋势 ,表现 出 明 


显 的 区 域 差 异 。 其 中 ,2003 一 2009 年 中 亚 冰川 表面 高 程 总 体 下 降 了 9.59+1.89 m;I 区 ( 即 西藏 和 青海 南部 ) 的 冰川 表 
面 高 程 下 降 了 6.51+2.9 m; I (BIAI 、 青 海北 部 和 甘肃 部 分 地 区 ) 下 降 了 7.87+5.03 m;II 区 ( 即 中 国境 外 ,中 亚 地 


区 的 部 分 国家 ) 下 降 了 9.81+5.1 m, 且 监测 到 2004 一 2005 年 冰川 表面 高 程 上 升 。 本 研究 方法 对 冰川 区 点 云 类 高 程 


脚 点 监测 具有 


定 的 通用 性 ,但 对 基准 DEM 的 依赖 度 较 高 。 


关键 词 : 冰川 表面 高 程 ; ICESat; 卫星 激光 测 高 ;时序 重 建 ; 高山 冰川 ; 中 亚 


在 全 球 变 暖 背景 下 ,气候 变化 加 剧 ,中 亚 冰 川 
变化 表现 出 明显 区 域 差异 。 青 藏 高 原 中 部 和 西北 
部 冰川 相对 稳定 ,而 高 原 周边 山区 的 冰川 物质 亏损 
严重 。 喀 喇 昆仑 地 区 冰川 异常 ,冰川 跃 动 现象 延 
伸 到 昆仑 西部 和 帕 米 尔 ,与 中 亚 其 他 地 区 的 负 质 量 
平衡 形成 鲜明 对 比 ”。 目 前 , 喀 喇 昆仑 冰川 质量 增 
加 暂缓 , 且 开 始 消融 ”。 冰 川 变化 形式 复杂 多 样 , 冰 
川 积 累 或 消融 必须 考虑 表面 高 程 变化 因素 。 了 解 
当前 和 未 来 冰川 质量 的 变化 ,以 此 预测 干旱 区 冰川 
融 水 径流 的 未 来 变化 情况 和 自然 灾害 的 发 生 , 提 前 
提出 防治 措施 ,对 于 淡水 资源 安全 与 经 济 可 持续 发 
展 等 有 重要 意义 “。 

目前 ,测定 冰川 变化 主要 有 花 杆 测量 法 .卫星 
定位 监测 方法 .激光 测 高 等 5”。 冰 川 区 坡度 大 、 冰 本 
徐 盖 和 冰 前 湖 发 育 特 征 使 冰川 变化 的 局 部 复杂 化 ， 
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增加 了 冰川 高 程 变 化 监测 的 难度 中。 随 着 遥感 技术 
的 兴起 ,由 于 其 多 时 相 大 范围 覆盖 的 特点 ,被 广泛 
应 用 于 测定 冰川 高 程 变化 。 

本 文采 用 2003 一 2009 年 的 ICESat 卫星 激光 测 
高 数据 ,构建 适合 中 亚 高 山区 冰川 表面 高 程 数据 的 
处 理 模 型 ,获取 中 亚 地 区 2003 年 各 脚 点 每 年 的 拟 合 
高 程 , 最 后 根据 拟 合 高 程 结 合 降水 .温度 以 及 地 形 
等 因素 分 析 研 究 区 冰川 表面 高 程 变 化 特征 及 驱动 
因素 。 


1 原理 与 方法 


1.1 点 云 去 品 及 其 精度 优化 方法 

ICESat 卫星 在 中 低 纬度 地 区 没有 精密 航天 顺 指 
向 控制 ,造成 单一 重复 轨道 在 中 低 纬度 地 区 并 不 完 
全 匹配 ,在 轨道 上 相隔 几 百 米 甚 至 几 公里 “。 本 文 
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通过 建立 去 噪 算法 和 多 项 式 模型 等 提高 ICESat 激 
光 测 高 数据 在 空间 和 时 间 上 的 一 致 性 。 

本 文 以 SRTM 高程 数据 作为 基准 面 ,采用 标准 
差 法 吻 除 由 于 受 云 层 等 因素 影响 的 ICESat 高程 异 
常 值 ,从 而 保证 实验 数据 的 可 靠 性 。 所 有 激光 脚 点 
与 基准 DEM 的 差 值 服从 正 态 分 布 时 , 设 激光 光束 脚 


点 与 基准 DEM 的 高 程 差 值 为 D=[d,d,, ds;,*……， 


dy_,, dy] ,其 平均 值 为 4 ,标准 差 为 e,d, 表示 第 i 个 
高 程 脚 点 ,N 为 激光 脚 点 的 总 个 数 。 则 激光 脚 点 与 
基准 DEM 的 高 程 差 值 落 在 (-3 e ,+3e) 的 概率 如 下 
式 所 示 。 
P(-3e<d,-d<+3e) ~99.7% (1) 
如 果 d, 服 从 正 态 分 布 ,在 3 售 标 准 差 的 方法 中 ， 
异常 值 就 可 定义 为 与 平均 值 偏差 超过 3 倍 标准 差 的 
值 ;如 果 d; 不 服从 正 态 分 布 ,可 以 采用 远离 平均 值 
的 多 倍 标准 差 来 描述 。 后 续 在 2.3.2 节 中 ,具体 说 明 
如 何 根 据 正 态 分 布 结果 ,采用 远离 1 倍 标准 差 ,来 易 


除 异 常 值 。 
1.2 大 范围 冰川 高 程 变化 时 序 重建 
D SRTM 高 程 为 自 变 量 ,采用 公式 (1) 噜 除 置 信 


度 低 的 高 程 脚 点 ,并 建立 区 域 高 程 拟 合 多 项 式 模 
型 ,确定 每 年 激光 雷达 卫星 ICESat 脚 点 高 程 与 
SRTM 高 程 之 间 的 关系 。 对 每 组 给 定 的 数据 拟 合 点 
(x> y) (其 中 ,i<m,m 为 脚 点 数量 ; x 为 中 脚 点 72 对 应 
的 重 采样 SRTM 高程; y 为 ICESat 激 光 光 束 脚 点 i 的 
高 程 ), 为 了 获取 x 与 y, 的 n(nsm) 次 多 项 式 函 数 关 
Z P(x))。 
P (x)= Sazi (2) 
式 中 :为 第 k(k<n) 多 项 式 ;a 为 拟 合 多 项 式 的 系 
数 。 当 拟 合 值 P,(xj) 与 y; 越 接近 ,表示 激光 脚 点 i 获 
取 的 高 程 值 越 精确 。 当 所 有 激光 脚 点 的 拟 合 高 程 
P (x) 与 ICESat 脚 点 高 程 y 差 值 的 平方 和 7 值 最 小 
时 ,由 系数 a, 构成 的 多 项 式 模型 也 就 越 精 确 。 
T=min (Ère) -可 = ¥ [Sau -| (3) 
当 拟 合 函 数 为 n 次 多 项 式 时 ,满足 式 (3) 的 
P e) 称 为 最 小 二 乘 拟 合 多 项 式 ; 当 n=1 时 称 为 线性 
WE, 为 系数 oo， …，aw (和 =0,1,…* ,nn) 的 多 
元 函数 。 上 述 问 题 即 为 求 T= T(ay,a,,°°*, a) 的 极 
值 问 题 。 对 公式 (3) 以 a 为 未 知 数 求 导数 为 0 时 
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(j=0, 1,.…, 
式 如 下 : 


m 


S [San -Bra eter a 
i=0 \k=0 i=0 


将 公式 (4) 按 照 关 于 a, (J=0,1……， n ER, 
用 和 矩阵 表示 为 式 (5) 的 法 方程 组 。 


m m 
met Şa om Yar] | Sy, 
is i=0 i 
m m m m 
$a Se Deel [Say 
i i i t| = iJi 
i=0 i=0 i=0 ea | iso (5) 


n), 得 多 元 函数 求 极 值 的 必要 条 件 公 


公式 (5) 的 系数 矩 阵 是 对 称 正定 矩阵 ,存在 唯 
一 解 。 当 n=3 时 ,求解 得 到 的 拟 合 高 程 P a) 与 基准 
高 程 x, 相关 系数 RS 0.95 , 选 n=3 得 到 区 域 的 拟 合 
高 程 , 即 为 激光 测 高 某 个 时 间 激 光 脚 点 i 的 拟 合 高 
程 值 与 基准 DEM 的 拟 合 多 项 式 如 下 : 

P (x) =a% +a,x; ta,x, + a, (6) 
式 中 : x, 为 基准 高 程 值 ; P,(%) 为 待 校正 高 程 脚 点 数 
据 ; a, , a, , a ,a6 分 别 为 拟 合 多 项 式 的 系数 。 


2 实验 


2.1 研究 区 概况 

中 亚 山 区 位 于 亚洲 中 部 海拔 在 3500 m 以 上 , 冰 
川 面积 约 1.3x10” km。 为 研究 中 亚 高 山 冰 川 区 内 部 
的 冰川 变化 差异 ,本 文 根 据 地 貌 、 地 形 等 生态 地 理 
区 域 系 统 将 研究 区 划分 为 3 个 区 域 ( 图 1)。I 区 主要 
包括 西藏 及 青海 南部 地 区 ,该 区 东南 部 冰川 大 部 分 
属于 海洋 性 冰川 ,内 部 冰川 大 多 属于 大 陆 性 冰川 ; 工 
区 主要 包括 新 疆 地 区 及 青海 北部 和 甘肃 北部 的 部 
分 区 域 , 该 区 冰川 以 大 陆 性 冰川 为 主 ,但 天 山 \ 昆 仑 
山东 段 以 及 喀 喇 昆仑 山北 坡 等 部 分 地 区 冰川 属于 
亚 大 陆 性 冰川 ; 王 区 在 中 国境 外 ,与 中 国 相 邻 ,地 热 
东 高 西 低 ,大 部 分 为 温带 大 陆 性 气候 。 
2.2 数据 与 预 处 理 

本 文 使 用 数据 主要 包含 以 下 3 部 分 : (1) ICESat 
数据 来 源 于 冰 、 云 和 陆地 测 高 卫星 ,卫星 于 2003 年 1 
月 13 日 成 功 发 射 , 是 首 颗 激光 测 高 卫星 ,搭载 地 学 
激光 测 高 系统 GLASN ;ICESat 数 据 产品 在 美国 国家 
冰雪 数据 中 心 (http://www.nside.org) 人 免费 获取 。 根 据 
研究 目的 ,本 文采 用 2003 一 2009 年 的 ICESat 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2016)1665 
号 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 中 亚 冰 川 分 区 概况 示意 图 


Fig. 1 Overview of glacier zoning in Central Asia 


GLAH14 数 据 , 为 便于 时 序 重建 ,避免 降雪 和 融化 等 
季节 间 冰 川 高 程 变化 的 影响 ,选择 每 年 的 同一 月 份 
数据 ;根据 统计 的 数据 量 , 选 择 数 据 量 最 多 日 最 完 
整 的 每 年 3 月 的 数据 进行 处 理 分 析 。(2) 本 文选 择 
能 够 免费 获取 的 90 mm 分 辩 率 包含 研究 区 范围 的 
SRTM3 数据 作为 参考 DEM。(3) 冰川 编目 数据 采用 
2017 年 7 月 发 布 的 RGI 6.0(Randolph Glacier Invento- 
ry 6.0) 数 据 库 ,RGI 6.0 以 2010 年 左右 的 夏季 无 云 
Landsat TM/ETM+ 为 主要 影像 来 源 , 用 于 获取 人 研究 
区 冰川 边界 。 
2.3 数据 处 理 

首先 从 ICESat 数 据 中 提取 脚 点 的 坐标 、 高 程 ， 
进行 坐标 转换 和 地 理 配 准 ; 再 噜 除 受 地 形 、 云 层 等 
因素 影响 产生 的 异常 值 ,本 文 根 据 正 态 分 布 检验 结 
果 选 定 异 常 值 并 将 其 剔除 。 利 用 去 除 异 常 值 后 的 
ICESat 坐标 数据 ,在 SRTM DEM 数据 中 重 采样 , 获 
取 ICESat 脚 点 坐标 对 应 的 SRTM DEM 高 程 ,建立 每 
年 的 多 项 式 拟 合 函数 ,以 解决 ICESat 数 据 重复 轨道 
不 完全 匹配 的 问题 ,并 计算 相关 系数 尼 , 作 为 评 佑 
拟 合 效果 参数 。 以 2003 年 的 SRTM DEM 高 程 数据 
作为 基准 ,获取 2003 年 每 个 脚 点 数据 每 年 的 拟 合 高 


程 ,实现 时 序 重 建 。 
2.3.1 数据 预 处 理 ”由 于 本 文 始终 以 SRTM DEM 为 
基准 进行 时 序 重 建 , 因 此 需 将 Topex/Poseidon 高 程 
转换 为 WGS84 高 程 ,如 下 式 。 

Hey = hgs - N-0.7 (7) 
式 中 : Hy 为 转换 后 点 云 WGS84 高 程 ; ha, 为 ICESat 
点 云 相对 于 Topex/Poseidon 参 考 椭 球 的 高 程 (m); M 
是 大 地 水 准 面 和 参考 椭 球 面 的 距离 (m) ;数值 0.7 的 
单位 为 m, 为 2 个 参考 椭 球 间 的 高 程 差异 ""。 本 文 
在 30 kmx30 km 尺度 范围 内 选取 最 高 点 作为 SRTM 
DEM 山顶 点 ,对 比 参 考 椭 球 为 WGCS84 高 分 辨 率 四 
维 地 球 影 像 筛选 均匀 分 布 可 靠 的 山顶 点 ,记录 其 作 
为 控制 点 的 150 个 山顶 点 的 经 纬度 坐标 ,对 SRTM 
DEM 进行 配 准 ,提高 不 同 坡度 . 坡 向 的 ICESat 脚 点 
RARE” 
2.3.2 DRA RA APRA ok) a Fe h OF 
究 区 内 ICESat 数 据 高 程 点 。 以 工区 的 2003 年 数据 
为 例 ,检验 ICESat 脚 点 高 程 是 否 符合 正 态 分 布 。 如 
表 1 所 示 , 显 著 性 为 0.00<0.05 ,所 以 该 区 域 高 程 值 不 
服从 正 态 分 布 , 因 此 , 选 定 远离 平均 值 1 倍 标 准 差 的 
值 为 异常 值 并 将 其 剔除 ,确保 数据 的 准确 性 。 
2.3.3 多 项 式 拟 合 函数 ”以 SRTM DEM 高 程 为 自 变 
量 , 与 ICESat 高 程 进行 多 项 式 拟 合 , 建 立 多 项 式 拟 
合 函 数 如 式 (6) ,并 计算 相关 系数 尼 。 如 图 2 所 示 ， 
通过 对 比 工 区 内 2003 年 点 云 数 据 去 除 异 常 值 前 后 
拟 合 效果 ,发 现 去 除 异 常 值 的 尼 更 接近 1, 拟 合 函 数 
的 拟 合 效果 更 好 ,更 符合 冰川 表面 高 程 的 实际 
情况 。 


3 ”结果 与 分 析 


3.1 多 项 式 拟 合 模 型 评估 

由 于 卫星 数据 的 处 理 过 程 带 有 保密 性 ,难以 给 
予 物理 模型 对 数据 源 进 行 校正 。 本 文通 过 统计 学 
的 方法 ,采用 确定 系数 RR 评估 多 项 式 拟 合 效果 , 表 
示 模 型 对 现实 数据 拟 合 的 程度 , R 的 计算 公式 


表 1 I 区 的 2003 年 ICESat 脚 点 高 程 正 态 分 布 检验 
Tab.1 Normal distribution test of ICESat foot elevation in area I in 2003 


Kolmogorov-Smirova 检验 


Shapiro-Wilk 检验 


统计 自由 度 显著 性 


统计 自由 度 显著 性 


0.041 1483 0.000 


0.983 1483 0.000 


注 :Kolmogorov-Smirova 表 示 标 化 后 的 数据 是 否 服从 理论 的 分 布 ;Shapiro-Wilk 表示 一 种 基于 相关 性 的 算法 ,显著 性 假设 检验 方法 。 
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图 2 工区 2003 年 数据 去 除 异 常 值 前 后 对 比 


Fig. 2 Comparison of data in zone I before and after removing outliers in 2003 


如 下 。 


2_ SSR _SST-SSE_， SSE 
人 =SST= SST al SST (8) 


式 中 :SSR 为 预测 数据 与 原始 数据 均值 之 差 的 平方 
和 ;SST 即 原始 数据 和 均值 之 差 的 平方 和 ;SSE 为 拟 
合 数据 和 原始 数据 对 应 点 的 误差 的 平方 和 ; R 的 取 
值 范围 为 [0,1]。 当 拟 合 数 据 和 原始 数据 对 应 点 的 
误差 越 小 时 ,SSE 越 小 , 拟 合 数据 越 可 靠 , 此 时 民 的 
取 值 越 接 近 1。 因 此 , 当 R 的 取 值 越 接近 1, 拟 合 效 
果 越 好 。 

由 表 2 所 示 ,整体 冰川 区 IT 区 JI 区 JI 区 的 R? 
都 大 于 0.98 ,说 明 整 体 拟 合 效果 很 好 。 冰 川 区 由 物 
质 平衡 线 分 为 消融 区 和 积累 区 ,在 物质 平衡 线 以 下 
为 消融 区 ,物质 平衡 线 以 上 为 积累 区 ,冰川 表面 高 
程 变化 可 总 体 归 为 这 2 种 模式 ,因此 利用 三 次 多 项 
式 更 能 符合 2 种 模式 的 冰川 表面 高 程 变化 现象 , 且 
符合 冰川 表面 高 程 自身 变化 规律 。 
3.2 冰川 表面 高 程 变化 分 析 

在 数据 处 理 的 过 程 中 由 于 系统 误差 会 产生 部 
分 异常 值 , 本 研究 采用 式 (9) ,根据 正 态 性 检验 结 
R ,将 数据 中 与 变化 平均 值 的 偏差 超过 2 倍 标准 差 
的 高 程 变化 平均 值 作 为 异常 值 进行 剔 除 。 各 区 不 
同时 间 段 内 冰川 表面 高 程 平均 变化 量 见 表 3。 


E +E +E? E; 
T= 1 = no ~ (9) 


式 中 :为 误差 = 高 程 变 化 值 -高 程 变化 均值 。 
由 图 3a 可 知 ,整体 区 域 即 整 个 中 亚 地 区 所 有 的 


5800 6200 6600 


SRTM 高 程 /m 


表 2 研究 区 2003 一 2009 年 3 月 多 项 式 拟 合 模 型 的 R 
Tab. 2 R’of the polynomial fitting model from 2003 to 
March 2009 in the study area 


年 份 整体 工区 工区 TK 

2003 0:9951 0.9866 0.9974 0.9916 
2004 0.9990 0.9995 0.9994 0.9959 
2005 0.9964 0.9924 0.9980 0.9963 
2006 0.9987 0.9979 0.9992 0.9975 
2007 0.9985 0.9994 0.9997 0.9932 
2008 0.9995 0.9988 0.9998 0.9989 
2009 0.9999 0.9999 0.9999 0.9998 


冰川 区 域 , 总 的 冰川 表面 高 程 下 降 9.59+1.89 m, ŒE 
2004 年 高 程 下 降 了 5.89+3.76 m, 次 年 回升 到 与 2003 
年 相当 的 高 度 , 之 后 逐年 缓慢 下 降 ; 图 3b 中 ,区 冰 
川 表面 高 程 在 2003 一 2004 年 间 下 降 5.78+2.13 m, 
2005 年 略 有 回升 ,之 后 高 程 逐 年 缓慢 下 降 , 变 化 斜 
率 为 -0.581 ma'; 图 3c 中 ,I 了 [区 其 冰川 表面 平均 高 
FETE 2003—2009 年间 下 降 了 7.87+5.03 m, 2005 年 
冰川 表面 高 程 出 现 大 幅 上 升 ,甚至 超过 2003 年 的 高 
程 值 ,但 整体 仍 处 于 下 降 趋 势 ,变化 斜率 为 -1.226 
ma'!; 图 3d 所 示 , 与 其 他 分 区 相 比 ,厂区 每 年 的 冰川 
表面 平均 高 程 变 化 波动 较 大 。 

如 图 4 所 示 , 整 个 研究 区 内 天 山 、 帕 米尔 高 原 和 
青藏 高 原 东 南部 地 区 的 冰川 高 程 下 降 剧 烈 ,青藏 高 
原 南 部 地 区 冰川 高 程 呈 轻 微 下 降 趋 势 。 中 部 的 昆 
仑 山西 段 地 区 部 分 冰川 高 程 呈现 升 高 状态 , 西 段 以 
及 昆仑 山东 部 的 祁连山 地 区 都 旦 下 降 趋 势 。 北 部 
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表 3 各 区 不 同时 间 段 内 冰川 表面 高 程 平均 变化 量 


Tab.3 Average variation of glacier surface elevation in different periods of time in each region 


时 段 整体 mm 工区 /mm 工区 /mm 亚 区 /mm 
2003—2004 4F —5.89+3.76 —5.78+2.13 —3.53+2.99 -4.82+3.31 
2004—2005 年 6.13+3.91 1.79+1.44 8.13+5.52 1.83+0.52 
2005 一 2006 年 -3.78+2.66 -1.12+1.01 -6.71+5.24 —2.02+1.68 
2006 一 2007 年 —3.56+2.79 -0.35+0.30 -4.61+3.31 -0.03+1.89 
2007 一 2008 年 —1.97+1.04 -0.63+0.60 0.43+0.32 —4+3.53 
2008 一 2009 年 —0.52+0.38 —0.42+0.29 —1.58+1.41 —0.77+0.64 
2003 一 2009 年 —9.59+1.89 -6.51+2.9 —7.87+5.03 —9.81+5.1 

(a) 整体 变化 : (b) 工区 变化 
5070 [ 一 国 一 平均 值 「 _ 一 国 一 平均 值 
[ — 趋势 线 S772. & = 趋势 线 
5067 上 | mm 5770.19 

r | 5066.28 y=-—1.4692x+8007.9 : 5770 i y= -0.58146x+6931.76 

aH 5064 上 a 5768 | 1 

tE | 5062.5 1E L 

R 5061 T R 

x F 5766 + 

P- L 5058.93 Bb L 

li + 「 5763.68 
5055 eS gg ES. ap 50 | 44 5762 Pog ,EE et ES ES ne EE SS | | | 
2003 2004 2005 2006 2007 2008 2009 2003 2004 2005 2006 2007 2008 2009 
年 份 年 份 
5025 r (©) 工区 变化 m- 平均 值 - (d 亚 区 变化 m- 平均 值 
l: — 趋势 线 4384 H T -一 趋势 线 
上 | | 4381.62 
5020 5018.12 y= —1.663x+7712.86 

£ j g 4380 上 4378.63 

a 5015 L T y= —1.226x+7470.3 ae | 4376.61 

1E 1E i 

g F 5011.41 g 4376 F 

BE 5010 S0023 ig | 4372.57 

F 一 一 4372 上 
5005 上 5006.8 | 
5005.66 4371.81 
到 l 1 l 1 L 1 | L | L 1 L L J 4368 l 1 l L 1 L 1 1 1 1 1 1 1 J 
2003 2004 2005 2006 2007 2008 2009 2003 2004 2005 2006 2007 2008 2009 
年 从 年 从 


图 3 各 个 研究 区 冰川 表面 高 程 变 化 时 序 重建 


Fig. 3 Time series reconstruction of glacier surface elevation change in each study area 


的 天 山地 区 冰川 高 程 变化 呈现 出 明显 的 区 域 差 异 ， ”区 .昆仑 山大 部 分 地 区 以 及 青藏 高 原 中 部 的 部 分 地 
东 天 山 . 北 天 山 和 西天 山南 部 地 区 冰川 高 程 下 降 明 区 冰川 呈 轻 微 增 厚 状态 ,但 本 文中 在 该 部 分 地 区 数 
显 ,西天 山北 部 冰川 高 程 下 降 速 率 相 对 较 慢 。 整 个 据 稀 少 , 无 法 确认 冰川 变化 情况 。 且 本 文中 颜色 机 
研究 区 内 除 昆仑 山 部 分 地 区 冰川 高 程 增 厚 外 ,其 他 格 表示 小 区 域 范围 内 激光 点 的 变化 率 , 因 此 ,与 
地 区 整体 旦 下 降 趋 势 , 只 是 下 降 速 率 存 在 地 区 差 。 “Treichler 等 中 研究 中 大 分 区 冰川 表面 高 程 变化 率 存 
异 。 与 2019 年 Treichler 等 ”的 研究 相 比 ,各 个 区 域 。 在 部 分 差异 。 

的 冰川 表面 高 程 变化 趋势 基本 符合 。 但 由 于 本 文 根据 2019 年 Treichler 等 2 对 MERRA-2 和 ERA- 
只 采用 了 每 年 3 月 的 数据 ,数据 量 小 ,分 布 稀 跑 ,部 。 ”Interim 气 候 再 分 析 数 据 的 处 理 ,显示 在 2003 一 2009 
分 区 域 缺少 数据 支持 且 表 面 高 程 变化 特征 表现 性 ” 年 间 , 气 温 呈 上 升 趋 势 , 除 昆仑 山 以 及 周 赎 少数 地 区 
差 , 例 如 在 Treichler 等 ”的 研究 中 ,中 天 山东 南部 地 外 ,亚洲 冰川 表面 高 程 整体 呈 下 降 趋势 。 青 藏 高 原 
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图 4 各 区 冰川 表面 变化 速率 分 布 


Fig.4 Distribution map of glacier surface change rate 


in each region 


EEKE ,整体 呈现 自 东 南 向 西北 递减 的 分 布 规 
律 中 ,东北 部 干旱 区 雨季 降水 量 也 呈 增 加 趋势 中 。 
图 3 中 ,在 2004 一 2005 年 间 各 区 的 冰川 表面 高 程 变 
化 剧烈 。 据 相关 研究 ”显示 ,2003 年 12 月 至 2004 
年 2 月 ,研究 区 冬 .春季 大 部 分 地 区 气温 较 常年 同期 
明显 偏 高 ,降水 偏 少 ,造成 冰川 消融 。2004 年 新 疆 
南部 等 地 夏季 高 温 日 数 普遍 达到 10~20 d ,局 部 地 区 
达到 40 qd"”。2005 年 西部 地 区 伴 有 降水 降雪 的 增 
加 ,造成 各 区 的 冰川 表面 高 程 在 2005 FAE ETH 
势 "9。 其 中 ,新疆 东 南部 和 西部 以 及 青藏 高 原 北部 
等 地 较 常 年 同期 偏 高 3 倍 ”。 

HJE E ,冰川 在 低 海 拔 坡度 小 于 30" 时 ,冰川 厚 
度 减 薄 随 着 坡度 的 减 小 而 加 剧 ,冰川 表面 高 程 降低 
较 快 ;中 高 海拔 地 区 ,温度 逐渐 降低 ,受降 水 的 影响 
较 大 , 易 导 致 冰 川 表面 高 程 增 高 ; 而 在 高 海拔 地 区 
平均 坡度 >35°, 易 发 生 冰 崩 、 雪 朋 等 现象 ,使 得 冰川 
物质 减少 中 。 结 合 地 形 因 素 ,降水 随 海拔 增高 而 增 
多 ,但 有 一 个 最 大 降水 量 高 度 ,超过 此 高 度 ,山地 降 
水 不 再 随 高 度 递增 , 随 气 候 干 湿 而 异 ”。I 区 中 , 青 
藏 高 原平 均 海拔 4000 m 以 上 ,海拔 高 度 对 气温 的 影 
响 已 超过 纬度 位 置 作用 。 因 此 ,青藏 高 原 除 南部 和 
东南 部 的 冰川 退缩 较 严 重 外 ,大 部 分 区 域 冰 川 高 程 
变化 相对 缓慢 下 。II[ 区 南部 有 昆仑 山 和 阿尔 金山 ， 
这 2 个 地 区 海拔 高 ,干旱 少雨 ,冰川 高 程 变化 相对 稳 
定 , 但 天 山西 北部 受降 水 量 影响 较 大 的 外 山脉 冰川 
融化 速度 较 快 。 祁 连 山地 表 平 均 温 度 在 海拔 高 的 
地 区 变 暖 趋势 明显 ,因此 ,祁连山 地 区 冰川 表面 高 
程 缓慢 下 降 辣 。 亚 区 中 西 帕 米尔 高 原 山 脉 交错 复 
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林 , 降 水 丰富 ,该 地 区 冰川 高 程 变 化 强烈 。 东 帕 米 
尔 高 原 东南 边缘 高 山 阻隔 了 来 自 印 度 洋 、 太 平 洋 的 
暧 湿 气 流 , 气 候 干 燥 , 冰 川 高 程 缓慢 下 降下。 就 下 
降 趋 势 而 言 ,研究 区 内 境外 区 域 冰川 表面 高 程 下 降 
趋势 最 快 ,为 -1.663 m.a ,新 疆 地 区 次 之 ,青藏 高 原 
地 区 最 为 缓慢 。 


4 结 论 

本 文 针 对 ICESat 卫星 测 高 数据 ,考虑 到 重复 轨 
道 在 中 低 纬 度 地 区 空间 匹配 性 弱 的 问题 ,选取 
SRTM DEM 为 高 程 基准 面 ,建立 冰川 区 点 云 去 噪 及 
其 精度 优化 算法 和 多 尺度 冰川 区 表面 高 程 时 间 序 
列 分 析 模 型 ,为 点 云 类 遥感 监测 冰川 表面 高 程 变化 
提供 参考 。 主 要 结论 如 下 : 

(1) 采用 正 态 分 布 显著 性 检验 ,进行 粗 差 剔 除 ， 
在 不 服从 正 态 分 布 情况 下 , 选 定 远离 线性 回归 1 倍 
标准 差 的 高 程 值 为 异常 值 ,能 够 有 效 噜 除 冰 川 区 IC- 
ESat 脚 点 高 程 中 的 异常 值 。 

(2) 利用 三 次 多 项 式 模 型 ,对 中 亚 山 区 整体 和 
分 区 域 的 不 同 尺度 冰川 表面 高 程 进行 拟 合 , 模 型 良 
都 在 0.98 以 上 ,表明 ICESat 数 据 和 SRTM 高 程 数 据 
的 三 次 多 项 式 关系 在 该 区 具有 较 强 的 普 适 性 。 

(3) 亚洲 高 山区 的 冰川 表面 高 程 2003 一 2009 
年 整体 下 降幅 度 在 10 m 以 内 ,表现 出 明显 的 区 域 差 
异 , 对 于 2004 年 前 后 表现 出 强烈 变化 ,由 于 环流 和 
季风 等 因素 的 影响 ,造成 的 温度 以 及 降水 异常 ,从 
而 导致 冰川 异常 波动 。 

本 研究 方法 对 冰川 区 点 云 类 高 程 脚 点 监测 有 具 
有 一 定 的 通用 性 ,但 会 使 点 云 数据 更 稀 玻 , 另 对 基 
准 DEM 的 依赖 度 较 高 。 后 续 , 研 究 将 进一步 加 入 其 
他 冰川 区 点 云 数 据 ,研究 时 间 序 列 建立 所 适用 的 空 
间 尺 度 , 建 立 更 灵活 多 变 的 区 域 冰川 表面 高 程 时 序 
监测 模型 ,对 影响 冰川 表面 高 程 变 化 的 地 形 、 环 境 
等 因子 进行 更 有 效 全 面 的 分 析 。 
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Abstract: The topography of the alpine glacier area in Central Asia is complex. Satellite laser altimetry 
technology can be used to monitor large-scale glacier surface elevation change. This study used ICESat laser altimetry 
data from 2003 to 2009, and shuttle radar topography mission (SRTM) elevation data and Randolph glacier 
inventory (RGI) 6.0 cataloging data for 2000 to monitor the surface elevation change of alpine glaciers in Central 
Asia. The glacial region of Central Asia is divided into three regions according to geographical form. First, the 
data were tested for normality to eliminate error. When a normal distribution was not observed, outlier values 
from the standard deviation of the linear regression were selected as abnormal. Point cloud denoising and its 
accuracy optimization algorithm and a multi-scale analysis model of temporal and spatial variation in surface 
elevation in the glacier area were established. A cubic polynomial model was used to fit the glacier surface 
elevation data for different scales across the whole area and sub regions of the high Asian mountains. A time 
series reconstruction of changes in alpine glacier surface elevation in Central Asia was carried out based on the 
data for each region from 2003. Temporal changes in glacier surface elevation from 2003 to 2009 were analyzed. 
The results showed that the average elevation for the whole region decreased, with clear regional differences. 
From 2003 to 2009, the surface elevation of high Asian glaciers decreased by 9.59 + 1.89 m. The glacier surface 
elevation in Areas I (i.e., Tibet and southern Qinghai), II (i.e., Xinjiang, Northern Qinghai, and some parts of 
Gansu), and III (i.e., outside China and some countries in Central Asia) decreased by 6.51 + 2.9 m, 7.87 + 5.03 m, 
and 9.81 + 5.1 m, respectively. The glacier surface elevation increased from 2004 to 2005. Area I showed the 
slowest decline in glacier elevation, followed by Area II; Area III showed the fastest decline. This research 
method has universal application potential for monitoring point cloud elevation in glacier areas. The model 
evaluation parameter (R’) was >0.98, indicating that the cubic polynomial relationship between ICESat data and 
SRTM elevation data showed strong universality in this area. However, this will make the point cloud data more 
sparse, and the model is highly dependent on the benchmark DEM. 

Keywords: glacier surface elevation; ICESat; satellite laser altimetry; temporal reconstruction; alpine glacier; 
Central Asia 


